Quasi-equilibrium during aging of the 2D Edwards- Anderson model 
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We test the quasi-equilibrium picture of the aging dynamics -strictly valid in the asymptotic 
dynamical regime of aging systems- in the pre-asymptotic aging regime of the two dimensional 
Edwards- Anderson spin glass model. We compare the fluctuation-dissipation characteristic for spin 
autocorrelation function and response with a corresponding one obtained for a suitably defined new 
correlation function and its conjugated response. In agreement with the quasi-equilibrium picture we 
find that after a short transient the two corresponding fluctuation-dissipation ratios (FDR) coincide 
at equal times. Moreover we show that, as it happens for the usual FDR, the new dynamic FDR at 
finite time coincides with the static one at finite size. 

PACS numbers: 05.70.Ln, 75.10.Nr, 75.40.Mg 



INTRODUCTION 

In recent times, following developments in spin glass 
mean field theory P, |3| , much emphasis has been put on 
the study of ofF-equiUbrium fluctuation-dissipation rela- 
tions during aging dynamics in glassy systems. These 
relations quantify the deviation of the ratio between cor- 
relation functions and conjugated responses from the one 
implied by the fluctuation dissipation theorem valid at 
equilibrium, and have been posed at the basis of a de- 
tailed thermodynamical and statistical description of the 
dynamics of glasses H, ^ |^ ■ Linear response theory al- 
lows to relate possible asymptotic violation of the fluctu- 
ation dissipation theorem to the failure of ergodicity at 
the level of the equilibrium measure [^. 

Given the correlation function C(t, t^) of a certain ob- 
servable A, and its conjugated response function xi^i'tw) 
describing the effect at time t of a fleld conjugated to 
A acting from time to time tw^ one can define the 
fluctuation-dissipation ratio (FDR) X{t,tw) from the re- 
lation: 



X{t,t^) = T 



dc{t,t^)/dt, 

dx{t,t^)/dt 



(1) 



This is just unity in equilibrium conditions while devi- 
ates from it off-equilibrium and in general depends on 
the observable quantity A at hand. In mean fleld spin 
glasses the FDR admits a non trivial limit in the aging 
regime, where the correlations assume a scaling invariant 
behaviour. Moreover, in the long time limit, taken after 
the thermodynamic limit, one can deflne the function 



x{q)= Yiu\ X{t,t^)\c(t.u,)= 



t,t 



(2) 



which can have a non trivial behavior. When this 
happens, x{q) turns out to have an important co- 
variance property under exchange of the given observ- 



able A chosen in the measure of correlation and re- 
sponse. If we have a correlation function CA(t,tw) 
corresponding to an observable A and a correlation 
function CB{t,tw) corresponding to an observable B, 
and we define an auxiliary limiting function qsiqA) = 
linit,t„^oo C'B(t,tu,)|cA(t,t™)=<?A' the following relation 
holds 



XBiqBiqA)) = XAiqA), 



(3) 



the meaning of which is that the functions XA{t, t^) and 
Xsit, t^) coincide asymptotically for equal times. More- 
over, in ref. 6] it has been shown that in a large class of 
finite dimensional systems with short range interaction 
x{q) defined in an out of equilibrium context is deeply 
related to the nature of the equilibrium free-energy land- 
scape. In fact, considering the overlap probability func- 
tion (OFF) P{q) describing the statistics at equilib- 
rium of the correlations of the observable A in two con- 
figurations chosen with Boltzmann weight, the linear re- 
sponse theory implies: 



xiq) = / dq'Piq'). 
Jo 



(4) 



This equality implies that either x{q) and P{q) are both 
non trivial or they are both trivial and could be taken 
as the starting point for an experimental measure of 
the equilibrium OFF from off-equilibrium dynamics. Of 
course it does not imply the existence of some short range 
system where it is verified non-trivially. Going through 
the derivation one realizes that Q expresses the com- 
mutation of the thermodynamic limit and the long time 
limit as far as certain susceptibilities are concerned. No- 
tice also that, thanks to equation Q expresses the 
fact that for two observables A and B, couples of states 
with identical qA also have identical qB, i.e. the func- 
tion qA[qB) defined in the dynamics describes the re- 
lation between different overlaps in equilibrium ergodic 
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components. This property has been shown to be deeply 
related to ultrametricity in '^l where it was found that 
the combination of relations Q and 10} implies ultra- 
metricity. 

The meaning of eq. Q has been clarified in Q where 
it has been discussed how x{q) can be related to the den- 
sity of metastable states -or quasi-states- with free-energy 
density slightly above the minimum, implying that quasi- 
states of equal free-energy are selected with equal proba- 
bility during the dynamical process. The identity |0J and 
the covariance property lO, allow to rationalize jSilSj 
interpretation of the ratios T/x{q) for different values of q 
as effective temperatures governing the exchanges of heat 
among slow modes evolving on waiting time dependent 
time scales 0. 

Effective temperatures dependent on q mean that while 
modes evolving on the same scale are in equilibrium with 
each other, heat exchanges between modes evolving on 
widely separated scales do not occur. It has been re- 
cently shown |lOl llH that in trap models where x{q) 
has a non trivial q dependence, but ultrametricity does 
not hold, different quantities define different FDR. A sit- 
uation where it would be difficult to identify the FDR's 
with effective temperatures. Conversely, J.-L. Barrat and 
L. Berthier studied Lennard- Jones models of glass- 
forming liquids where an FDR constant in q seem to de- 
scribe the off-equilibrium dynamics, and found that den- 
sity fluctuations at different wave vectors define the same 
FDR. 

Numerical simulations of 3 and 4 dimensional spin 
glass Edwards-Anderson models, comparing extrapola- 
tions of the OFF from finite size systems and extrap- 
olations of the FDR from finite time indicate the non- 
triviality -and consistently the identity- of both functions 
[l3j |. This has been taken as an evidence in favor of a 
"replica symmetry breaking scenario" (RSB) for finite 
dimensional spin glasses. 

These extrapolations however have been questioned in 
a series of papers showing that the OFF in systems with- 
out RSB, i.e where the OFF is trivial in the thermody- 
namic limit, can be plagued by severe finite size effects 
such that for relatively small systems it appears simi- 
lar to what one expect for systems with RSB In 
such conditions RSB could be wrongly inferred from ex- 
trapolations of the finite size OFF P{q, L) of systems of 
too small sizes L, while the true P{q) is a trivial single 
(5-function. In the same way one could doubt that off- 
equilibrium dynamic simulation times are too short to 
reliably extrapolate the asymptotic FDR from the finite 
tw, and that the true asymptotic one is just a single flat 
step as in domain growth problems . 

On the experimental side it is clear that many systems 
with slow aging dynamics are found in pre-asymptotic 
regimes. A common phenomenon is the one of inter- 
rupted aging, found e.g. in where a slow dynamical 
regime similar to usual aging eventually crosses over to 



equilibrium behavior. In addition, even in three dimen- 
sional spin glasses, the paradigmatic systems where aging 
could persist indefinitely, one sees that many quantities 
are far from their final values. In particular the exper- 
iments of Ocio and Herisson [T^l where the first exper- 
imental determination of the FDR in spin glasses was 
achieved, show FD curves that strongly depend on the 
waiting time, signifying that the dynamics is still in some 
pre-asymptotic regime. In such conditions it is of great 
interest to inquire if the concepts valid for aging in the 
asymptotic regime can be adapted to get an adequate 
picture of the dynamics on much shorter time regimes. 

In this context, one can hypothesize that the identity 
between static and dynamic FDR found in 3 and 4 dimen- 
sion is due to the fact that at a given time tw there is a 
slowly growing length £,{tu,) over which the system has ef- 
fectively equilibrated, and X{q^ tw) would approximately 
respect the relation with P{q, L = ^{tw))- Such an ex- 
tension would suggest the approximate validity at finite 
time of a quasi-equilibrium picture of the aging dynamics 
in which quasi-states with equal free-energy are selected 
with equal probability, and the static-dynamic equiva- 
lence would just reflect the properties of the equilibrium 
landscape of finite size system. The hypothesis is rather 
suggestive as it would provide a framework to interpret 
aging properties in an appropriate time scale even for 
systems which display interrupted aging. Here, while in 
a certain time window slow evolution and approximated 
scaling laws for correlations and/or susceptibilities are 
observed, the final asymptotic state is ergodic. 

To test this extension L. Berthier and A. Barrat 
pit studied the 2D Edwards-Anderson (2DEA) model, 
which, on one hand displays strong aging effects at finite 
times, and a non trivial OFF for finite size, on the other 
it is known to finally reach a paramagnetic state at all 
finite temperatures. In that work it was found that in- 
deed it exists a correspondence L ^ t^ such that the 
relation holds. More recently, Berthier also studied 
the three and four dimensional case in a pre-asymptotic 
regime obtaining similar results [lit ■ 

In the light of the previous considerations about the 
link among effective temperatures and time scale sepa- 
ration, these finding appear rather surprising. Here, no 
time scale separation is possible, slow modes have to ex- 
change heat in order to eventually equilibrate. In or- 
der to save the picture, one can of course hypothesize 
that this exchange occur "adiabatically" , in such a way 
modes evolving at the same rate appear able to equili- 
brate at their effective temperature with faster or slower 
modes before exchange heat. If this consideration applies 
FDR corresponding to different quantities should appear 
approximately equal one to another. In order to test 
this hypothesis we consider as in the 2D Edwards- 
Anderson (2DEA) model where as mentioned aging is 
interrupted after a finite relaxation time. In our analysis 
we define some suitable correlation and response func- 
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tion not obviously related to the usual spin autocorrela- 
tion and its associated response, and compare the FDR 
for both couples of functions. In addition, in two dimen- 
sions we test the equivalence between the static and the 
dynamic FDR for the new quantities. 

Our results are then compared with analogous mea- 
sures in the Viana-Bray diluted spin glass, where Q is 
known to hold non-trivially. 

The remaining of the paper is organized as follows. 
In the next section, we introduce the relevant quantities, 
then we present and discuss the results of the simulations 
and finally the conclusions are outlined. 



DEFINITION OF THE OBSERVABLES. 

The model we will consider consists in a pair of spin 
glass systems with independent random coupling and 
identical number of spins coupled through random in- 
teractions. Before explicitly introduce the model let's 
say a few words to motivate this choice, keeping in 
mind that our task will be to compare FDR's corre- 
sponding to different correlation-response couples. In 
spin models, defined in terms of an exchange Hamil- 
tonian H = ^ll^j JijSiSj the natural and most com- 
monly used choice to probe dynamical correlations is 
the spin autocorrelation function at different times: 
C{t,tw) — N~^Yli{^ii^)Si{iw)) the corresponding "zero 
field cooled" susceptibility with respect to small local 
i.i.d. Gaussian fields hi with variance /i^, introduced in 
the systems at time tw and kept on at later times reads. 



Wh^^i {hiSi{t)) where the over-line denotes 



the average over the field. A second common choice is the 
"energy correlation function" also known as "link over- 
lap" CE{t,t^) = N-^Y.^,JJ^J{S^{t)SJ{t^)), and the as- 



sociated response: XE{t,tw) = -m^^iiJijhjSiit)). In 
mean field, for Gaussian long range Jij 's one can show 
that in the thermodynamic limit, choosing the variance 
of the Jy's to be equal to \/N one has for all times and 
with no assumption about the dynamics: 



C{t,t^ 
C{t, 



(5) 
(6) 



distributed, but independent quenched disorder and cou- 
pled by a random field Ri . The Hamiltonian of this com- 
pound system is defined by: 



2J 



where J/^ and Jfj represent the quenched disorder in 
copies 1 and 2 respectively and are quenched variables re- 
specting the lattice topology and otherwise taken as i.i.d. 
from a Gaussian distribution with mean and variance 
1. The variables i?j which couple spins with identical 
label in the two copies have been chosen randomly with 
values Ri — ±K. 

The dynamical spin autocorrelation function now 
reads: 

C{t,t^) = (2iV)-i J2 {S|m{t^.) + Sf{t)Sfit^)), 

i 

(8) 

and the corresponding response 

x{t, = ^ E (^'^'(^^ + w) (9) 

As second couple of correlation response pair we consider 
the spin cross-correlation function Q , 



Ccrossit,t^,) = (2iV)-i J2 i(SHt)S?{t^) + Sf{t)Sl{t^))R, 

(10) 



and 



ossit, t^) = — i-- ((^'"Slit) + hlSf{t))R,) 



2NhlK 

(11) 

Where the (...) indicates an average over the initial con- 
ditions and over the disorder. We will speak about 
direct correlation and response respectively for (jS)) and 
© and cross correlation and response for 11U|I and lll|l . 
An explicit formula for Ccrossit, tw) and Rcross{t, <«,) = 

— ^2<£j£^dLhd as functional of C(t,tj^) and R{t,t^) — 

— ^^^^^hA can be given for small K using linear response 
theory: 



so that, automatically, for all times, the FDR defined 
with these quantities coincide with the one defined with 
the usual correlations and response. Analogously at equi- 
librium, one finds that the relation qsiQ) = 9^ holds and 
2qPE{qE{q)) = P{q) independently of the ultrametric 
nature of the organization of the states. 

Then, in order to test the quasi-equilibrium picture 
one needs to compare overlaps non-trivially related one 
to the other. Consider, therefore, two copies of spin glass 
systems, with identical number of spins, and identically 



Ccrossit, t^) = K^I3[ f ™ dsCit, s)Rit^,s) (12) 
Jo 

+ ds C{t^,s)R{t,s) ] 
Jo 

Rcross{t,t^) = ds R{t,s)R{s,t^), 

which shows that if the cross FDR coincide with the di- 
rect one, it is for non trivial reasons. 
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RESULTS AND DISCUSSION 

We studied the cross-quantities in two different sys- 
tems of Ising spins with random quenched disorder, the 
Edward Anderson model, a bi-dimensional square lattice 
of spins of size N — L x L and the fixed connectivity ver- 
sion [i^l of the Viana-Bray model where the spins 
are on a random lattice with fixed connectivity c = 10 
and size N. For both models, we considered two copies 
with identical number of spins and independent quenched 
disorder coupled by a random field Ri = ±K as in the 
previous section with K — 1/2. For this large value of K 
we are out of the linear response regime that allowed us 
to derive the explicit form of the cross-quantities as func- 
tion of the usual ones, but even if the relations (|12|l does 
not hold, there is no reason to believe that the relation 
between X^ross and X becomes trivial. 

In order to have as a reference results for a system 
where the picture sketched in the introduction certainly 
holds we present first the data of quick simulations of the 
Viana-Bray model. 

In a first test, we compared the FD plots in dynamic 
simulations of aging experiments with the static ones ob- 
tained through the parallel tempering technique. Our 
results are summarized in figure ^ We can observe that, 
on one hand, the static characteristics has small finite size 
dependence, on the other, the dynamical curves little de- 
pendence on the waiting times. As expected the dynamic 
curves coincide with the static ones for the direct func- 
tions. For the cross functions they also coincide provide 
the static curves are shifted vertically, this is normal as it 
should be noted that, the maximum value of qcross in the 
equilibrium OFF of finite systems is unity, while dynam- 
ically Ccross{t,tw) is monotonically decreasing from the 
value Ccrossitw, iui) < 1 for t > t^,. Therefore, one should 
subtract a constant to the second integral of P{qcross)- 
S (qcross) — Xcrossiq')dq' to Compare it with the 

the dynamic function. 

We then tested to what extent equation (0) is valid 
when finite systems in statics are compared to systems 
evolved for finite aging times in dynamics. From the 
statics we get the functions P{q) and P (qcross) where q 
is defined as the usual (direct) overlap between two in- 
dependent replicas S = iS^Sf) and S' = iSl',Sf') : 

q{S,S') = ^E^S|Sl' + SfSf, while qcross{S,S') = 

mT,^R^iSlSf' + SfS}'). Then using ^ we derive 
the equilibrium quantities x{q) and x{qcross) and com- 
pare them with the results obtained from the dynamic 
simulations. 

Although we did not try to measure the joined prob- 
abihty P{q, qcross), we could extract a function 

Hcross (q) 

as implied by the relation ij^J and compare with the one 
directly obtained from the dynamics. The results can be 
seen in figure |21 where one can see that the static and the 
dynamic curves approach each other for values of qcross 



smaller than Ccrossitw,tw)- This is what one should ex- 
pect because , as discussed above, in dynamics this is 
the largest value of Ccross{t,tw), and tends to its hmit 
from below for tw — > oo. Conversely in statics, for finite 
systems, the probability distribution always extends to 
values of qcross larger than the maximum value for an 
infinite system. 



o.s 




0.2 0.4 0.6 o.s 1 



C(r,tJ,Ccr„Jt,iJ 

FIG. 1: Fluctuation dissipation plot for the VB model. A 
vertical shift to superimpose the equilibrium curve to the dy- 
namic one. N = 169, 196, 256 and 324 and = 10, 10^, 10^ 
and 10*, r = 2.18 




0.2 0.4 0.6 0.8 1 

C(t,tJ,q 

FIG. 2: Parametric curves of the cross correlation as a 
function of the direct one. The lines are the equilibrium 
curves = 144 and 196 and the points the dynamic ones 
t„ = 10, 10^10^ and 10*. We see that the curves approach 
each other for values of Qcrosa smaller than Ccroas{tw,tw), 
which seems to have reached its ^ oo limit. 

We then pass to the study of the two dimensional sys- 
tem. We studied the aging dynamics of the 2DEA model 
with unitary Gaussian couplings at T = 0.43, where no 
sign of thermalization can be observed in the correlation 
function up to waiting times as high as — 10^. In our 
simulations we used = 0.02 and checked to be in the 
linear response regime using the value = 0.01. On the 
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equilibrium side, using the parallel tempering technique 
[iz] we were able to calculate with good precision, the 
spin glass order parameter P(g), as well as the function 

Crossl- 
in figure 121 we present a fluctuation-dissipation plot 
for the model. As observed by L. Bcrthier and A. Bar- 
rat [l^ . one can superimpose the finite time curves of 
the direct functions with the second integral of the equi- 
librium OFF of suitable size systems at the same tem- 
perature. The points represent the data obtained by the 
dynamic simulation, while the lines those obtained study- 
ing the static of the model using the parallel tempering 
technique. We present for clarity only plots at two dif- 
ferent waiting times, — 10^ and t^, = 10'^, for the 
dynamic simulations and two lattice sizes, L ~ 8 and 
L = 10 for the static ones. The upper curves correspond 
to the usual functions (similar curves were already pre- 
sented in ref. 18]) , while the lower ones, reflects the 
cross functions. Unfortunately, it turns out that in order 
to superimpose the static and dynamic curves, a vertical 
shift in S{q) is not enough and an horizontal shift, should 
also be performed. 

0.4 I , , , F , , , 1 
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FIG. 3: Fluctuation dissipation plot for the EA model at tem- 
perature T = 0.43 and waiting times = 10^, 10'^ in compar- 
ison with equilibrium functions for systems of size L = 8, 10. 
A vertical and an horizontal shift are necessary to superim- 
pose the equilibrium and the dynamic curves. 

In order to understand this point we notice that differ- 
ently to what happens for the direct function for which by 
construction C{t,t) = 1, the value of Ccross{t,t) evolves 
in time. In such conditions the definition of the FDR 
in terms of the simple correlation is not necessarily the 
most appropriate. In fact, the study of running away 
systems (e.g. brownian motion or particles in non con- 
fining random potentials shows that a better def- 
inition is obtained considering the following combination 
of the correlation function Brr qs.ti t, t^) = \ [Ccross{i, t) + 
Ccross{tw,tw) - 2Ccross{t,tui)]UM- This obviously is not 
the only combination which could be used; e.g. in |ll| 
it was suggested the use of Fcross{t,tw) = Ccrossit,t) — 



Ccross{t,tw), however, we preferred to present the data 
using Bcross , which is symmetric in the configurations at 
the two times and admits the interpretation of a distance 
between the configurations at times t^^ and t. We checked 
all the results below using Fcross instead of Bcross which 
turn out to be equivalent within our numerical accuracy. 

0.3 I , , , , , , , , 1 




FIG. 4: Fluctuation dissipation plot for the 2D EA model 
using B and Bcross as abscissas. The agreement of the dy- 
namical characteristic and the static one for the cross quantity 
is comparable to the direct one. In this case no shift of the 
curves is needed. 
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B(l,tJ , b 

FIG. 5: Comparison between the relation between Bcross and 
B as obtained directly from the dynamics and supposing eq. 
^ in statics. Solid lines: the static data for L = 8, 10. Points: 
the dynamical data are = 10^10^ The sizes of the static 
data are L = 8, 10. 

The static analogous of the function Bcross{t,tw) 
is, given two configurations S and S', the quan- 
tity bcross{^T^') = 2 ^) 9^035(8', S') — 

2<Zcross(S, S')]. In order to compare the dynamical FDR 
with the static one we should be cautious of the symme- 
try of the Hamiltonian under contemporary reversal of all 
the spins. As discussed in the proper static probability 
distribution to compare with the dynamics is not the full 
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distribution but the distribution modulo the symmetry 
of the Hamihonian. For the function P{q), symmetric 
under q ^ —q one just needs to consider the positive q 
part of the function multiplied by 2 (for normalization). 
In order to eliminate this symmetry in the distribution 
P{bcross) in an analogous way, we can just consider the 
histogram of the hcross corresponding to configurations S 
and S' such that qcross{^, S') is positive. The results are 
presented in figure0]which shows that the agreement be- 
tween the static and dynamic curves for the cross quanti- 
ties is comparable to the direct ones, for the same lattice 
sizes. 

Next we compare the relation between Bcrossit,tw) 
and B(t,tw) obtained directly in dynamics and relating 
the values with equal x in statics as explained above. 
This is shown in figure where we show that for the 
times and lengths considered there is a good correspon- 
dence between statics and dynamics. 

Finally, given the good quality of our susceptibility 
data, we could take the derivatives of the x versus B 
characteristics so to compare directly at equal times the 
cross FDR Xcross{t, tw) with the direct one X{t, t^). The 
comparison is shown in figure El which shows that even 
for waiting times as short as tw = 10^ the two quantities 
are very close to each other. Preliminary results indicate 
that this is also the case in three and four dimension. 




t 



FIG. 6: Direct comparison of the FDR for the cross and 
direct quantities in the 2D EA model. The temperature is 
T = 0.43 and the waiting times are, from bottom to top 
tn, ~ 10^, 10^, 10^. After a short transient both quantities do 
coincide. The glass width is L=512 and averages are made 
upon 190 samples. 



CONCLUSIONS 

In this paper we have proposed the use of cross cor- 
relation functions and response in disordered systems to 
probe the existence of effective temperatures during ag- 
ing of disordered systems. We have compared the be- 



haviour of a mean-field model with the one of a paramag- 
net, where aging is a transient behavior. We find that the 
at equal time the FDR for direct and cross quantities co- 
incide within numerical error after a short transient. The 
correspondence among dynamical FDR at finite time and 
static one for finite size is confirmed as far as the cross 
quantities are concerned. These two findings support on 
one hand the idea that aging dynamics can be described 
in terms of effective temperatures, on the other that these 
temperatures are related to the density of states of finite 
systems on a scale L(tw)- Preliminary results in three and 
four dimension indicate that the same kind of behaviour 
is found in these systems for time scales much shorter 
than the ones needed to reach an asymptotic state. 
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